Field induced spin ice like orders in spin liquid Tb2Ti207 
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We have studied the field induced magnetic structures in the spin Uquid Tb2Ti207, in a wide temperature 
(0.3<r<270 K) and field {0<H<1 T) range, by single crystal neutron diff'raction with H // [110] axis. We 
combined unpolarized neutron data with polarized ones, analyzed within the local susceptibility model. A 
ferromagnetic-like structure with k = propagation vector is induced, whose local order at low field and low 
temperature is akin to spin ice. The four Tb ions separate in a and j3 chains having diff'erent values of the 
magnetic moments, which is quantitatively explained by taking the crystal field anisotropy into account. Above 
2 T and below 2 K, an antiferromagnetic-like structure with k = (0,0,1) is induced besides the k = structure. 
It shows a reentrant behavior and extends over a finite length scale. It occurs together with a broadening of the 
nuclear peaks, which suggests a field induced distortion and magnetostriction eff'ect. 

PACS numbers: 71.27.+a, 75.25.+Z, 61.05.fg 



Geometrical frustration now attracts considerable interest, 
as being a possible tool for tuning several physical proper- 
ties concomitantly. The pyrochlore lattice of corner sharing 
tetrahedra offers the best model of this frustration in three di- 
mensions, leading to short range magnetic orders such as spin 
liquids, spin ices and spin glasses without chemical disorder. 
R2Ti207 pyrochlores where R=Dy or Ho are model spin ices^ , 
with a ground state entropy akin to that of real ice^. Ho or 
Dy magnetic moments interact via ferromagnetic first neigh- 
bor interactions and are constrained to point along their local 
<111> Ising axes. The local spin ice structure, with 6-fold 
degeneracy, consists of two spins pointing in and two out of 
each tetrahedron. 

The spin ice degeneracy is lifted by a magnetic field H, 
leading to original magnetic transitions, which depend on 
the orientation of H with respect to the local < 1 1 1 > trigo- 
nal axes^. With H // [111], magnetization plateau and liquid- 
gas transition occur"^-^-^, providing the first evidence of mag- 
netic monopoles^. Applying H // [100] yields the first mag- 
netic example^ of the Kasteleyn transition known in polymers. 
With H // [110], the lattice divides into a and chains^' 
with diflTerent angles between H and the local Ising axes. 

In Tb2Ti207, the crystal field (CF) anisotropy is weaker 
than in canonical spin ices, and ferromagnetic (F) and antifer- 
romagnetic (AF) first neighbor interactions nearly compen- 
sate. This widely studied spin liquid^^'^^, where short range 
correlated moments fluctuate down to 50 mK, was also called 
a "quantum spin ice"^"^. Its magnetic ground state (GS) is very 
sensitive to perturbations. Lattice expansion induced by sub- 
stituting Sn for Ti yields a long range "ordered spin ice" in 
Tb2Sn207^^. An AF long range order (LRO) is induced under 
applied stress and/or magnetic field^^'^'^'^^. 

Up to now, the complex field induced magnetic struc- 
tures were not studied precisely in Tb2Ti207. We have com- 
bined high accuracy neutron diffraction experiments with both 
polarized and unpolarized neutrons, in a wide temperature 
(0.3<r<270 K) and field (0<H<1 T) range, to investigate 
them in detail. The polarized neutron technique, newly used 
for such compounds, allows very small moments to be mea- 
sured with great accuracy, which is crucial for low fields and 



high temperatures. It yields microscopic information on the 
magnitude and orientation of the field induced ordered mo- 
ments. Our data analysis is based on the local susceptibil- 
ity approach We show that it is valid in fields up to 1 T, 
and temperatures 5<r<270 K. It takes advantages of the high 
symmetry of the pyrochlore lattice which allows us to refine 
these complex magnetic structures with only two parameters, 
for any orientation of the magnetic field. 

For all temperatures and fields H // [1 10], we observe an F- 
like structure with propagation vector k = 0, characterized by 
magnetic Bragg peaks of the face centered cubic lattice. We 
show that it involves a and yS-chains as in model spin ices. But 
in Tb2Ti207, these chains have very diflferent moment values. 
This original eflTect is described quantitatively using the CF 
parameters derived in Ref,20. It provides a nice example of 
diflTerent moments induced by the field on the same crystallo- 
graphic sites of a homogeneous and highly symmetric lattice. 

Below 2 K and for fields above 2 T, we also observe a 
second family of Bragg peaks, belonging to the simple cu- 
bic lattice They are indexed in the cubic unit cell of Fd3m 
space group with a propagation vector k = (0,0,1). We studied 
this AF-like structure down to 0.3 K in fields up to 7 T, us- 
ing unpolarized neutrons. We discuss its observation together 
with the F-like structure and propose an explanation for its ori- 
gin. The whole analysis provides a microscopic description of 
the field induced magnetic structures in Tb2Ti207. 

The neutron diflTraction studies were performed on the 
diflTractometer Super-6T2 at the Orphee reactor of the Labora- 
toire Leon Brillouin^^. The field was applied close to a <1 10> 
axis, with a misorientation of 5°. We used unpolarized neu- 
trons of wavelength An = 0.90 A, collecting 200 reflections 
for each (H, T) set, for 10 temperatures between 0.3 K and 
50 K and 7 fields between and 7 T. We also used polarized 
neutrons of incident wavelength An = 0.84 A and polariza- 
tion Po = 0.91. We collected flipping ratios at 200 to 400 
Bragg peaks for each (//, T) set, for 8 temperatures between 5 
and 270 K, in a field of 1 T. The programs FULLPROF^^ and 
CHILSQ^^ were used to refine the magnetic intensities and 
flipping ratios, respectively. 
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FIG. 1: (Color online)Tb2Ti207 : Local anisotropic susceptibility el- 
lipsoids ;^/yr, measured at 10 K (a) and 270 K (b). Ellipsoids were 
scaled by temperature to compensate the Curie behavior of the sus- 
ceptibility, (c) susceptibility components;^// and;^^ versus T. The 
lines are calculations using the CF parameters of Tb2Ti207 (see 
text); d: Measured versus calculated flipping ratio at 10 K; the re- 
finement carried out with H=l T // [110] on 423 flipping ratios gave 
=0.939(15) yt/^/T and;t'i2=0.535(ll) i^bIT with the goodness of fit 

;r'=2.92. 



A single crystal of Tb2Ti207 was grown by the floating- 
zone technique. Its crystal structure was refined at 50 K 
within the Fd3m space group (434 reflections measured in 
zero-field), in agreement with powder data^^. 

In the pyrochlore lattice, the 16 Tb atoms of the cubic 
unit cell occupy a single site (16 J) with local symmetry 3m. 
They can be subdivided into 4 diflTerent groups according to 
the directions of their local <111> anisotropy axes. For H 
// [110] the Fdd2 group must be used in the magnetic mo- 
ment refinements^. It is the highest symmetry subgroup of 
the Fd3m space group which leaves the magnetization invari- 
ant. In this subgroup the I6d sites split into two diflTerent sub- 
sets corresponding to the Of-chains (where H makes an angle 
of 35.3° with the local <111> easy axis) and to the yS-chains 
(where it is perpendicular), respectively. 

In the local susceptibility approach the atomic site sus- 
ceptibility tensor accounts for the linear paramagnetic re- 
sponse of the moments to an applied field of arbitrary direc- 
tion. The induced magnetic moment at the 16 J site of the 
unit cell writes: =1^^H, where is a tensor of rank 3x3, 
whose components depend on the symmetry of the atomic 
site. The magnetic moments are in general not parallel to the 
field, and the magnetic structure is not necessarily collinear, as 
reflected by the non diagonal terms in the susceptibility ten- 
sor. In the cubic axes, the symmetry constraints on the suscep- 
tibility tensor for a magnetic atom at the 16 J site in the Fd?>m 
group imply: xn = X22 = X33', Xn = Xu = X23' So only 
two independent parameters need to be determined regardless 
of the field direction. The susceptibility components;^// and 
X±, respectively parallel or perpendicular to the local <111> 



easy axis, measure the anisotropy on a given Tb ion. They are 
given by the relations: x//=Xn-^^Xi2 ^"^^^Xi-^Xw-Xn- 

We determined the susceptibility components by perform- 
ing a least square fit of the site susceptibility model to the 
flipping ratios measured for H=l T // [110]. The refinement 
at 10 K (Fig. Id) shows the excellent quality of the fit. The 
magnetic ellipsoids determined from the site susceptibility pa- 
rameters are shown in Fig. la and lb at 10 K and 270 K. For 
a given ellipsoid, the radius vector gives the magnitude and 
direction of the magnetic moment induced by a field rotat- 
ing in space. The elongation of the ellipsoid increases with 
decreasing T, reflecting the increasing CF anisotropy, as the 
magnetic moments evolve from Heisenberg to Ising behavior, 
being more and more constrained to align along their local 
<111> axes. At 10 K the two components;^// and;^^ (Fig. 
Ic) along the easy and hard local axes diflfer by a factor 5. 
They still diflTer up to 200 K, but become practically equal at 
270 K. The lines in Fig[T]: were calculated using the CF inter- 
action in Tb2Ti207^^, with a magnetic field of 1 T applied ei- 
ther parallel or perpendicular to the local < 1 1 1 > axis. Good 
agreement with the data requires the AF exchange to be taken 
into account: using a self-consistent calculation, valid down to 
-10 K, we obtain a molecular field constant X - -0.35 T/yUg , 
close to the value derived in Ref,20. 

Once the site susceptibility parameters are known, one can 
easily calculate the magnitude and the direction of the mo- 
ments induced on each Tb site by a field applied in an arbitrary 
direction. This approach is valid for fields of 1 T and tempera- 
tures 5<r<270 K, a range where the susceptibility parameters 
obey a Curie- Weiss law (Fig. Ic). Below ~5 K, the local sus- 
ceptibility approach is no longer valid and we used unpolar- 
ized neutron diflTraction to determine the magnetic structure. 
We refined the unpolarized neutron data assuming 4 diflferent 
moments (in magnitude and orientation) on the 4 independent 
Tb sites of the cubic unit cell. The values plotted in Fig. 2 for 
H = 1 T // [1 10] combine the analysis from polarized neutrons, 
for 5<r<270 K, and unpolarized neutrons for 0.3<r<10 K. 
The two sets of data show a perfect overlap. 

For H=l T // [110], the k = structure (Fig. 2a) at 0.3 K 
clearly shows a local spin ice order, with two moments point- 
ing in and two out of a given tetrahedron. All moments keep 
close to their <111> easy axis. The four tetrahedra of the 
unit cell are the same, as in Tb2Sn207 ordered spin ice^^. 
But in striking contrast with standard spin ice structures, here 
the Tb moments have very diflTerent magnitudes (Fig. 2b). 
The a-moments reach 5.7(3) \ib at 0.3 K, close to the 5.7(2) 
yUfi value deduced from the crystal field study whereas the 
/^-moments only reach 1-2 yUg. The angular dependencies (not 
shown here), show that cK-moments keep close to their local 
< 1 1 1 > easy axis up to about 3 K, then reorient smoothly along 
H with increasing T . The yS-moments remain perpendicular to 
H up to about 3 K, then reorient abruptly along H. 

The presence of two types of moments {fix and ^2) with 
diflTerent values can be attributed to the misorientation of H 
with respect to the [110] axis (5°). At 1 T and 0.3 K, our 
refinements yield yS-moment values of 0.8(2) and 1.3(2) yUg. In 
contrast, the two types of a-moments keep almost the same 
values. 



3 




E 4 



Q 2 



a chains >■ ^ 


\ 

1 X 

1 10 100^ 


(3 chains 







10 100 
T(K) 



FIG. 2: (Color online)Tb2Ti207 : (a) field induced k = F-like struc- 
ture at 0.3 K (H=1T//[110]). Balls and arrows represent Tb ions and 
moments, respectively, (b) temperature dependence of the Tb mo- 
ments from unpolarized (• for of-chains and o for /^-chains) and po- 
larized (■ for Qf-chains and □ for /^-chains) neutron data. The two 
moments arise from the 5° field misorientation. The thick red lines 
are a crystal field calculation with a field H=\ T applied at 35.3° 
(resp. 90°) from the local [111] axis for a (resp. /?) sites. Inset: mag- 
netization M versus temperature: (a) from the neutron data, (A) from 
bulk measurements of Ref .20.. Thin solid lines are guides to the eye. 



The induced a and yS-moments vary with temperature as 
predicted by the crystal field level scheme of Tb2Ti207^^, in 
the whole T range 0.3-270 K (Fig. 2b). In the calculation, the 
field of 1 T was taken exactly along [110] axis, so one needs 
to average the measured yS-moments to compare with the cal- 
culated ones. The agreement is quite good. It shows that the 
non-zero yS-moments found in the k = structure, not seen in 
model spin ices^^ result from the weaker Tb anisotropy and 
not from the field misorientation. A more sophisticated model 
should involve both applied and internal fields, and their pre- 
cise orientations with respect to the crystal axes. We note 
that diflTerent fii and j32 moments are predicted by a CF cal- 
culation taking the field misorientation into account. At low 
temperature, the strongly Ising character should render the JS- 
moments much more sensitive to the field orientation than the 
Qf-moments (Fig. 15 of ref. ^2^, as observed. Knowing the a 
and yS-moments, one can also calculate the average magneti- 
zation M in one tetrahedron, (Fig. 2b inset), which perfectly 
agrees with the bulk magnetization data^^. 

When H increases, the moment values in the k = structure 
increase (Fig. 4a), and their angles with H decrease, showing 
that they reorient along the field. The F of-moments quickly 
saturate, whereas the F yS-moments increase more slowly with 
H, reaching 3.2(3) jub and 6.0(3) jjb fory^i andyS2, respectively, 
at 7 T and 0.3 K. 

For H > 2 T and below 2 K, new magnetic peaks appear, 
much weaker than the previous ones, indexed in the space 
group Fd3m with k=(0,0,l). We refined this new family of 
peaks as a second magnetic structure, also involving the whole 
sample, and analyzed independently. The k=(0,0,l) value 
means that in the cubic cell, two tetrahedra have the same mo- 
ment orientations and the other two are reversed, so that this 
AF-like structure has no net magnetization. In the k=(0,0,l) 
structure, only the AF yS-moments are significant, the AF a- 
moments are negligible (Fig. 3a and 3b). At 0.3 K, up to 
H=l T, the AF y^^-moments remain oriented at 10(5)° from 
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FIG. 3: (Color online)Tb2Ti207 : (a) field induced k = (0,0,1) AF- 
like structure at 0.3 K, (H=7 T // [110]). (b) field dependence of the 
moments in the AF-like structure; the a moments (not plotted) are 
below 0.2 (c) (H, Taf) transition line (•) our measurements, (A) 
from Ref 18. (d) AF correlation length versus H\ in inset, the (112) 
peak for several fields. Solid lines are fits with a Lorentzian peak 
shape, and dashed line is the resolution peak shape. 
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FIG. 4: (Color online)Tb2Ti207 : (a) temperature dependence of the 
Qf-moments in the F-like structure in high fields {H >2 T). The F a- 
moments increase with H and decrease with T below ~2 K when the 
AF-like structure settles in. (b) Field dependence of the full width 
half maximum (FWHM) of the (440) peak with negligible magnetic 
contribution. The value at H-Q is the resolution limit. 



their local <1 1 1> axis, and 90(5)° from the field. 

The k=(0,0,l) structure is stabilized in a certain (//, T) 
range, determined by plotting the peak integrated intensities 
versus T or H. This yields the transition line Taf(H) plot- 
ted in Fig. 3c. It agrees with a previous determination^^ for 
0.3<r<1.5 K, but deviates above. The k=(0,0,l) structure 
has a finite correlation length ^af, as shown by the Lorentzian 
linewidth of the magnetic peaks (Fig. 3d), ^af starts in- 
creasing with H, saturates at about 70 A around 4 T, then de- 
creases. Its variation reflects that of the AF y(5-moments (Fig. 
3 a), suggesting a reentrant behavior. This AF-like structure 
should disappear at very high fields, where all Tb moments 
are aligned. 

We now summarize the field induced GS in Tb2Ti207. Be- 
low 2 T, the GS consists of a LRO k=0 structure, with field 
induced moments both on a and yS-chains. The moments val- 
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ues are well explained by CF calculations. The Tb anisotropy, 
much weaker than in model spin ices, yields non-zero mo- 
ments on the /^-chains. Above 2 T, the GS is a superposition 
of two modes with k=0 and (0,0,1), the latter ordering with 
a finite length scale. The k= (0,0,1) structure occurs together 
with a decrease of the F a-moments in the k=0 structure (Fig. 
4a) and the onset of well defined spin waves It is also con- 
nected with a field broadening of the nuclear Bragg peaks 
(Fig. 4b), suggesting a lattice distortion. All these features 
reflect a symmetry breaking with respect to the quantum spin 
ice state stable at if =0. 

The k= (0,0,1) structure induced by H// [110] strongly re- 
sembles that induced by a stress along the same direction^^. 
This strongly suggests that this structure is induced by mag- 
netostriction efl'ects. For //=0, a small distortion was ob- 
served in Tb2Ti207 below 20 K, likely precursor of a Jahn- 
Teller transition^^. In applied field, taking the bulk modulus^^ 
(Bo=156 GPa) and volume magnetostriction^^ of Tb2Ti207 at 
4 K, we estimate that a field of 7 T is equivalent to a pres- 
sure of 0.03 GPa, well below the stress of 0.2 GPa that in- 
duces the AF order However this estimation based on the 
isotropic volume is only a lower estimate for the stress. In 
Tb2Ti207, magnetostriction coefficients // and _L H have very 
high values (up to 5. 10""^), of opposite signs which nearly 



compensate in the isotropic volume magnetostriction. The 
strong anistropic magnetostriction arises from the specific Tb 
crystal field. The sample length // H expands whereas that 
_L H contracts^^. At a microscopic scale, one could speculate 
that bonds in a-chains expand, weakening the near neighbor 
exchange interaction along these chains, whereas bonds in 
yS-chains contract, reinforcing This could explain the AF 
field induced periodicity along the yS-chains. 

In conclusion, we determined the field induced magnetic 
structures in Tb2Ti207 by combining polarized and unpolar- 
ized neutron diff'raction, in a wide temperature and field range 
(H // [1 10]). The low field GS consists in a k=0 (F-like) struc- 
ture, with the local structure of a spin ice, but with Tb mo- 
ments of very diff'erent magnitudes. With increasing field, the 
Tb moments reorient from the spin ice easy axes to the field 
axis. Above 2 T, a k=(0,0,l) (AF-like) structure is also sta- 
bilized. It exists in a limited (7, H) range, has a finite length 
scale and shows reentrant behavior. It is attributed to mag- 
netostriction. As temperature increases, the susceptibility el- 
lipsoids show a progressive change from Ising to Heisenberg 
behavior. This change is demonstrated using the local sus- 
ceptibility approach, and well explained by the crystal field 
anisotropy of the Tb ion. 
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